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Abstract 

We study pattern formation associated with the polarization degree of free- 
dom of the electric field amplitude in a mean field model describing a nonlin- 
ear Kerr medium close to a two-photon resonance, placed inside a ring cavity 
with flat mirrors and driven by a coherent x-polarized plane-wave field. In 
the self-focusing case, for negative detunings the pattern arises naturally from 
a codimension two bifurcation. For a critical value of the field intensity there 
are two wave numbers that become unstable simultaneously, corresponding to 
two Turing-like instabilities. Considered alone, one of the instabilities would 
originate a linearly polarized hexagonal pattern whereas the other instabil- 
ity is of pure vectorial origin and would give rise to an elliptically polarized 
stripe pattern. We show that the competition between the two wavenum- 
bers can originate different structures, being the detuning a natural selection 
parameter. 
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I. INTRODUCTION 



Spatiotemporal patterns in the transverse direction of an optical field have now been 
widely studied theoretically and experimentally |0J. Studies of optical pattern formation 
share a number of aspects and techniques with general investigations of pattern formation 
in other physical systems |J, but they also have specific features such as the role of light 
diffraction and the vectorial degree of freedom associated with the polarization of the light 
electric field amplitude. A prototype simple model which has been very useful for the under- 
standing of pattern formation in nonlinear optical cavities is a mean field model describing 
a Kerr medium in a cavity with flat mirrors and driven by a coherent plane- wave field [|3],|J . 
This model was extended to take into account the polarization degrees of freedom in |5|-0. 
Some of the basic polarization mechanisms of pattern formation in alkali vapors or other 
non-linear materials can be understood in terms of this simple model despite the fact that 
the model is too simple to give a complete description of alkali vapors. Furthermore, the 
relative simplicity of the model in |J makes it worthwhile to study it in depth as a general 
prototype model for the basic understanding of vectorial patterns. A first study was under- 
taken in || for the case of linearly polarized driving field and the positive cavity detuning. A 
more detailed study in which the case of elliptically polarized driving field is also considered 
has been presented recently in reference J7J. 

Cavityless Kerr media can be used as an optical phase conjugation mirror. An ideal 
phase conjugate mirror should generate an output field such that the amplitude of the 
field, its propagation vector and its polarization unit vector are the complex conjugates of 
the corresponding magnitudes in the input field. The first two properties can be achieved 
easily using, for example, four- wave mixing processes 0. Usually the third property can be 
obtained only when the pump waves are circularly polarized and counterrotating. In order to 
find a system that verifies the three properties (vector phase conjugation) for an arbitrarily 
polarized input field it is required to make use of the special tensor properties of two-photon 
atomic transitions in degenerate four- wave mixing processes ||[10|. More precisely, vector 
phase conjugation can be achieved if the two-levels coupled by the two-photon transition 
have equal angular momenta J with J = or J = 1/2. In this situation the X1122 component 
of the susceptibility tensor vanishes. Intuitively, as AJ = 0, the atom do not change its 
angular momentum either by absorbing two pump photons or by emitting a probe and a 
conjugate photon so that the conjugate photon must be emitted with angular momentum 
equal and opposite to that of the probe photon. A detailed calculation || shows that indeed 
this is true for the cases indicated before. We should stress that the polarization properties 
of two-photon-resonant degenerate four-wave mixing processes are different from those of 
most other degenerate four-wave mixing processes. In the two-photon case the underlying 
physical mechanism is scattering of the probe field from a spatially uniform temporally 
varying coherence induced by the two pump waves whereas in the other cases is scattering 
from a spatially varying refractive-index distribution induced by the interference between 
the pump and the prove beams. Experimentally, vector phase conjugation was first observed 
using the 3Si/ 2 — > 6S1/2 two-photon transition in sodium vapor |TT . 



Here we will show that for a ring cavity filled with a nonlinear Kerr medium close to a 
two-photon-resonance (so that X1122 = 0) and illuminated with linearly polarized input and 
with negative detuning a new interesting situation appears: A codimension two bifurcation 
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in which two stationary Turing-like instabilities occur simultaneously. The first instability, 
if the other were not present, would originate a hexagonal pattern which is polarized linearly 
in the same direction as the input field. The second instability is of pure vectorial origin and 
if the first instability were not present it would give rise to an elliptically polarized stripe 
pattern. Here we study the interplay between the two instabilities. The codimension two 
bifurcation appears here in a natural way associated to the two-photon-resonant four-wave 
mixing nonlinearities rather that as the result of the fine tuning of two system parameters 
as is usually the case. The intensity of the pump field is the single control parameter to be 
tuned to change the distance to the instabilities. Furthermore, the system still have another 
easily accessible control parameter, the detuning, which allows the system to form different 
patterns while remaining at the same distance to the codimension two instability threshold. 
In particular we show how the detuning can be used as a tuning parameter to select the 
pattern. 

The outline of this paper is as follows: In Sec. [TT] we describe the model we are considering, 
its spatially homogeneous solution and the stability analysis of this solution. In Sec. |I 1 1 
using a weakly nonlinear analysis we derive the evolution equation for the patterns arising 
from the interacting instabilities. From these equations, the selected patterns for different 
values of the detuning are analyzed in Sec. |TV|. In Sec. [V] we describe the results from 
numerical simulations of the model and finally in Sec. [VI| we give some concluding remarks. 



II. DESCRIPTION OF THE MODEL, REFERENCE STEADY STATES AND 

STABILITY ANALYSIS 

The system we consider is a ring cavity filled with an isotropic Kerr medium. The cavity 
is driven by an external x polarized input field. The situation in which the polarization 
degree of freedom of the electromagnetic field is frozen was first considered by L. A. Lugiato 
and R. Lefever ||£|]. Geddes et al || generalized the model of [|J to allow for the vector 
nature of the field. Their description of this system is given by a pair of coupled equations 
for the evolution of the two circularly polarized components of the field envelope E + and 
E_, defined by 

E ± = -j=(E x ±iE v ). 
For an isotropic medium, the equations are 

^± _ n i ;„a\ tp , „-„v72i7 ,77 , • r a\ tp 12 , / a , dm m 121 



dt 



-(1 + i V e)E ± + iaV 2 E ± + E + ir][A\E ± \ 2 + {A + B)\E T \ 2 }E±, (1) 



where Eq represents the components of the input field (the right and left circularly po- 
larized components are equal since we consider x linearly polarized input), rj = +1 (— 1) 
indicates self-focusing (self-defocusing), 9 is the cavity detuning, a represents the strength of 
diffraction and V 2 is the transverse Laplacian. The parameters A and B are related to the 
nonlinear susceptibility tensor components in the following way: A = 6x1122 and B = 6x1221 
||. Also, for an isotropic medium we have A + B/2 = 1. As discussed in the introduction, 
here we are considering two-photon transitions between levels with equal angular momenta 
J where J = or J = 1/2, so that X1122 = (A = and B = 2). Also as we consider the 
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self- focusing situation, in what follows we take r\ = +1. The intensity of the input field is 
J = 2\E \ 2 . 

The steady state homogeneous solutions of Eq. ([I]) are reference states from which 
transverse patterns emerge as they become unstable. There is a symmetric (E s+ = E s _ = 
Eg) and two asymmetric (E s+ ^ E S J) homogeneous solutions 0. The symmetric solution 
corresponds to linearly polarized output light, while the asymmetric solutions to elliptically 
polarized output. Increasing the input field, the asymmetric solutions appear only for values 
of Iq larger than the instability threshold for pattern formation J7|, so here we will only 
consider the symmetric solution \ T2 \, 

I /2 = I s (l + (2I s -9) 2 ), (2) 

which gives an implicit formula for I s = \E S \ 2 . As it is well known, Eq. (0) implies bistability 
for 9 > y/3. However here we will always consider negative detunings which are far away 
from the bistable regime. An example of the symmetric solution for linearly polarized input 
is given in the inset of Fig. |l| 

Basic features of the stability of the steady state homogeneous symmetric solution can 
be analyzed by considering the evolution equations for perturbations ip± defined by 

E ± = E s (1 + 4j ± ). (3) 

From Eqs. (U) and (f|) we find 

dtip± = -[1 + i(0 - 2I S ) - zaV 2 ty± + tI s 2(^ T + r T + I^tI 2 )(1 + V>±) • (4) 

It is convenient to make a change of variables to the following basis 

E = K a 2 , a 3 , a,f = + ^+ + ^+ ~ ^-), ^+ ~ ^-)) T , (5) 

where T stands for transpose. In this basis, which emphasizes the role of symmetric = 
and antisymmetric = —if>-) modes, Eq. (|]) may be written as: 

<9 4 £ = L£ + iV 2 (£|£)+iV3(£|£|£), (6) 

where the linear matrix L is a matrix with 2x2 blocks in which the symmetric (ci, 02) and 
antisymmetric (03, (T4) modes are decoupled. 

As a consequence, the linear instabilities lead to the growth of either a symmetric or an 
antisymmetric mode. In Fourier space we have: 

f -1 (6 - 21. + ak 2 ) \ 



and 



-1 (6-2I s + ak 2 ) 

+ 2I s + ak 2 ) -1 



(9) 
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where k = \k\. 

Instability occurs if at least one of the eigenvalues A of L\ and L 2 has a positive real 
part. In Fourier space, these eigenvalues are solutions of the characteristic equations: 



(Ai + 1) 2 
(A 2 + l) 2 



+ ak 2 -6I S )(9 + ak 2 - 2I S 
+ ak 2 + 2I s )(9 + ak 2 - 2I S 







(10) 



For 9 < 0, we have 



A 1± = -1 ± ^4/2 - cfi(k 2 - k*y 
A 2 ± = -1 ± v / 4/ s 2 -a 2 (A; 2 -A:2)2, 



(11) 



where ak 2 = 4J S + \ 9\ and ak 2 = \9\. For both eigenvalues, the instability occurs at I s = 1/2. 



The non linearities in (||) include quadratic iV 2 (£|£) and cubic terms A^ 3 (E|S|S): 



N 2 {E\E) 



( 2a 3 a^ — 2a 1 a 2 \ 

3(Ti(Ti + O 20" 2 — O 3CT 3 + (J4(T4 
2cr 2 (J3 — 2<Ti<T4 
— 20i0 3 — 20"20"4 / 



(12) 



and 



AT 3 (£|£|£) 



2 



^ 2criCr3Cr 4 — cr 2 (o"lO"l + 0"20"2 — 0"3°"3 + CT4d"4) \ 

a 1 ((j 1 (T 1 + o 2 o 2 + cr 3 cr 3 — cT4cr 4 ) — 2o 2 o 3 o 4 
2ai<j 2 a 3 — a A (aiai — o 2 o 2 + a 3 a 3 + cr 4 cr 4 ) 
\a 3 (—aiai + a 2 o"2 + <J 3 a 3 + 0-40-4) — 2o!0 2 o 4 J 



(13) 



The structure of these terms also gives some general information on the nature of the insta- 
bilities. In particular, if the quadratic non linearity iV 2 (£|£) does not vanish, one expects 
the formation of hexagonal patterns instead of stripes. As explained in Ref. when the 
symmetric mode becomes unstable an hexagonal pattern is expected whereas when the an- 
tisymmetric mode becomes unstable there are no relevant quadratic couplings so a stripe 
pattern is formed. 

In Fig. ^ we plot marginal stability curves for 9 = — 1 as a function of ak 2 . The shape 
of the marginal stability curves is, in fact, the same for any value of the detuning 9. This is 
because the eigenvalues Aj given by Eq. flTTp depend on ak 2 — \9\, so a change in the value 
of 9 is equivalent to a displacement of the origin of ak 2 by the same amount. The origin 
moves to the right if the detuning 9 is increased. 

The instability region I comes from the eigenvalue Ai + so the critical modes are sym- 
metric and of zero frequency. A subcritical hexagonal pattern is expected via a transcritical 
bifurcation. If this were the only instability, it would correspond to the case discussed in [|J , 
in which the polarization degree of freedom is not taken into account. This instability leads 
to an x-polarized pattern while the ^-polarized component of the field continues to be zero. 

The instability region II comes from the eigenvalue A 2+ so the critical modes are anti- 
symmetric and of zero frequency. A stripe pattern is expected M. Given the antisymmetric 
nature of the unstable mode, the x-polarized component of the field is stable and remains 
almost homogeneous, while the stripe pattern appears in the y-polarized component, which 
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has zero value below the instability. Overall, the electric field displays an elliptically polar- 
ized spatial structure. We remark that such an instability is of pure vectorial nature with 
no analogue when the polarization degree of freedom is frozen. 

In the case considered here, 9 < 0, starting from the linearly polarized homogeneous 
solution, as the input field is increased, the system crosses the two instability thresholds 
simultaneously. This is a codimension two bifurcation involving two sets of stationary modes. 
The critical modes associated to region I are symmetric and have a critical wave number k s , 
while the critical modes associated to region II are antisymmetric and have a critical wave 
number k a . The ratio k a /k s can be changed easily varying the value of 9. 



III. WEAKLY NONLINEAR ANALYSIS FOR INTERACTING TURING 

INSTABILITIES 

The eigenmodes of the linear evolution matrix L are, in Fourier space (S + (k), 0, 0, 0) T , 
(0, 0,0) T , (0,0,i+(£),0) T and (0, 0, 0, i_(£)) T , with 

S ± (k) = ± (3 s (k)a 2 (k), A ± (k) = a 3 (k) =F (3 a (k)a 4 (k), (14) 

where U(k) = J r [U](k) denotes the Fourier transform of U(r). Furthermore 

UP S - a 2 (k 2 - k 2 ) 2 



where index j stands for s or a and /3j(kj) = 1. The critical modes correspond to the 
eigenvalues Xi + (k s ) and A 2 +(fc a ). Note also that 

S + + S_ S+-S- + . A+-A- 
ai = ^— ' a2 = —2pT> a3 = —^> a4 = 2&~~- (16) 

After diagonalization of the linear evolution matrix, the dynamics @ may be rewritten, 
in Fourier space, as 

d t S(k) = A(k)S(k) + N 2 (S\S)\- k + N 3 (S\S\S)\ n , (17) 

where S = (S + , S^, A + , AJ) T and A is a diagonal matrix with diagonal elements 
{Ai + , Ai_, A 2 +, A 2 _}. N 2 \k an d N 3 \^ are the Fourier transforms of the nonlinear terms of 
equation ([]), where the <Tj have been replaced by the corresponding linear combinations of 
S + , A + and A_. Slightly above threshold, this dynamics may be reduced to the dynam- 
ics of the critical modes amplitudes only, through the adiabatic elimination of the stable 



non-critical modes. This procedure is now standard fLo] , and we will only sketch here the 
main steps of its application to model (|]), and derive evolution equations for the critical 
modes, up to cubic nonlinearities. 

Let us write the amplitudes of the critical modes as S{k) = S + (k)5(\k\ — k s ) and A{k) = 
A + (/c)5(|A;| — k a ). Their dynamics writes 
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S(k) = (21. - l)S(k) + (N 2 \- kik=ks ) s + (N 3 \ k , k=ke )s , 
A(k) = (2I S - l)A(k) + (N 2 \ % , k=ka )A + (N 3 \ k , k=ka )A. 



(18) 



The terms in the quadratic nonlinearities N 2 are convolutions of products of critical and 
non-critical modes. The terms involving non-critical modes only (St, At, 5_A_) may be 
neglected, since they contribute, through the adiabatic elimination process, to nonlinearities 
of quartic or higher order. Furthermore, the structure of the cubic nonlinearities of Eq. @ 
is such that the pure critical mode contribution to N 3 vanishes. Hence N 3 will be neglected. 
We are thus dealing here with a non-generic case, where cubic nonlinearities are generated 
solely by the adiabatic elimination of stable modes from quadratic terms. The nonlinearities 
{N 2 \ k k=ks )s an d (N 2 \^ k=ks )A can be written at the lowest order in a non-critical modes as 



{N, 



k 
2 



dk' 



S(k - k')S{k') - A{k - k')A(k') 



{N, 



1\k,k=kJ A 



+1 
I 



+2S(k - k')S-(k')8(\k'\ - k 8 ) - 2A(k - k')A„(k')5(\k'\ - k a ) 
S(k)(S + (0) + £_(0)) + S(-k)(S+(2k) + S42k)) 

dk' \s{k-k')A{k') 

+S{k - k')A^(k')5(\k'\ - k a ) + A(k - k')S4k')S(\k'\ - k s ) 
+I s A(k) [s + (0) + £_(0) 



(19) 



The non-critical modes S±(0), S±(2k 8 ), S-(k a ) and A_(k a ) present in the Eqs. fll9|) may 
be expressed as an expansion in powers of critical ones using the adiabatic elimination 
procedure. One has at the leading order 



5±(0) = 


Is f ' 

2Ai,±(0) J 1 


S±(2k s ) = 


L (1 


2A li± (2A; s ) 1 


S-{k s ) = 


/. f 


2\ 1 ,-{k s ) J 


A-(k a ) = 


- L 1 





dk 



A(k)A(-k) + (1 T 2/3 a (0))S(k)S(-k) 



lT2p s (2k s ))S(k s )S(k s 



dk' 



3S(k s - k')S(k') + A(k s - k')A(k') 



dk'S(k a - k')A(k') 



(20) 



The substitution of Eqs. ( ^Of) in Eqs. (|I9|) leads to the following asymptotic dynamics for 
the critical modes, valid close to the instability threshold 



S(k) = (2I S - l)S(k) + J dk' \v S(k - k')S(k') + Vl A(k - k')A(k') 
dk' J dk"u(k, k', k")S(k - k')S(k' - k")S(k") 
dk' ! dk"w(k, k', k")S(k-k')A(k' - k")A(k") 
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A(k) = (21, - l)A(k) + v 2 J dk'S(k-k')A(k') 

- J dk' J dk"w(k, k', k")A(k - k')A(k' - k")A{k") 

- j dk' J dk"u(k, k', k")A(k - k')S(k' - k")S(k") , (21) 

where v = —V\ — y = 4f , an d 

u(k, k', k") = u^dk'l - k s ) + u 2 5(k') + u 3 8(k')S(k" - k) 

w{k, k', k") = w^dk'l - k a ) + w 2 5(k') + w 3 5{\k-k' - k"\ - k a ) 

u{k, k', k") = Ui5(\k'\ — k s ) + m 2 5(A; / ) — w 3 5(\k — k'\ — k a ) 

w{k, k', k") = w 1 5{\k'\ - k s ) + w 2 5(k') , (22) 

3/^ 9 34-2 1 0\ /^13-|-6|$| 1^ 1^ 

with u x = u 2 = 4(2^, u 3 = -f (2^p, wi = w 3 = and W 2 = (^T^. 

It has to be noted that the u terms come from quadratic resonances between critical and 
non critical symmetric modes, while the w terms come from quadratic resonances between 
symmetric and antisymmetric modes. As in the case of isolated Turinglike instabilities, 
there is no quadratic resonance between critical antisymmetric modes only. As a result, 
pattern formation is expected to strongly depend on the existence of quadratic resonances 
between symmetric and antisymmetric modes. Hence, for the sake of simplicity, we will 
consider separately the case with quadratic resonance between symmetric modes only, and 
the case with quadratic couplings between antisymmetric and symmetric modes, which is 
more intricate. In the latter case, quadratic couplings are such that k\ = k 2 + k 3 , with 
\ki\ = k s and \k 2>3 \ = k a . Since coscfi = j^r, where <fi is the angle between symmetric and 
antisymmetric vectors, quadratic resonance arise only for \9\ > ^f. 



IV. PATTERN SELECTION AND STABILITY 

In this section we study the various patterns that may appear as asymptotic solutions 
of eq. (|2"TD, and their stability. Each of these patterns is built on an arbitrary number of 
critical modes pairs, and it is the nature of their nonlinear couplings that determines their 
stability. In the following, we label the modes contributing to the formation of a pattern 
S{ki) = Si, A{ki) = Ai. We also use the notation Si for the complex conjugate of Si. 



A.k a <!f or \9\ < ^ 

In this case, there is no quadratic coupling between symmetric and antisymmetric modes, 
and there is no contribution coming form the terms with the coefficients V\, v 2 , W\ and ^3. 
Let us then consider separately amplitude equations for each type of modes. 



S 



1. Antisymmetric modes. 



For patterns built on antisymmetric modes only, the amplitude equations for an arbitrary 
number m of pairs of modes are: 

m 

A, = (21, - l)A t - 2w 2 A t I A, | 2 . (23) 

j'=i 

Hence, a pattern built on m pairs of wavevectors is marginally stable versus a m + 1 pair 
of wavevectors. So, at this level of analysis, any pattern with an arbitrary number of 
wavevectors is possible, including patterns of the form AJ (k a r), where Jo is zeroth order 

Bessel function and A = \J ^ 2I 2 ~^ ■ 

The linear growth rate of the evolution of symmetric modes in the presence of such 
patterns is zero, so antisymmetric patterns are marginally stable versus symmetric mode 
patterns. It would be necessary to go to higher orders in the amplitude equations to complete 
the pattern selection analysis in this case. 



2. Symmetric modes. 



On the other hand, for patterns built on symmetric modes only, the amplitude equations 
for a triplet of such modes are: 

Si = (2I S - l)Si + 2v S i+1 S i . 1 - (2u 2 + u 3 )\S t \ 2 S t - 2{u x + u 2 )S i ^ \Sj\ 2 . (24) 

At the instability threshold (I s = |), hexagonal pattern appear via a subcritical bifurcation. 

Increasing the value of I s , stripes may also become stable for 21 s — 1 > 9(2+^)2 • There is a 
small region of bistability of stripes and hexagons and for larger values of I s only the stripes 
remain stable. 

Hexagonal patterns are stable versus antisymmetric modes for I s < | |l + ^^q^p 
Hence, sufficiently close to threshold, one may expect hexagonal patterns. On increasing I s , 
such patterns should become unstable. 



3. Mixed modes. 



The amplitude equations for mixed structure formed by triplets of symmetric modes and 
an arbitrary number of antisymmetric modes are: 

Si = (21 s - l)Si + 2v S i+ iSi-i 
-S 



A k = (21 a - l)A k - A k 



(2u 2 + m 3 )|^| 2 + 2(ui + u 2 ) \Sj\ 2 + 2w 2 \ A i 

j=i 1=1 

3 m 
3=1 1=1 



(25) 
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However, these equations do not admit non trivial steady states. As a result, in these 
conditions, hexagonal or striped patterns of symmetric modes, and patterns built on an 
arbitrary number of antisymmetric modes may be simultaneously stable. 



B. < k a < k s or \9\ > ^ 

In this case, quadratic resonances may occur between symmetric and antisymmetric 
modes, and one may now expect contributions coming from the V\, v 2 , Wi and W3 in the 
amplitude equations, which have to be modified accordingly. Let us then consider the 
different types of patterns that may arise in this case, and which are built on modes belonging 
to the following set of critical modes (up to an arbitrary phase angle) (cf. figure ||). 



1. Symmetric modes. 

Striped and hexagonal symmetric modes patterns are now always unstable versus anti- 
symmetric modes. This is due to the positive nonlinear renormalization induced by these 
modes in the evolution equation for their resonantly coupled antisymmetric modes (this 
renormalization is associated to the terms with — w% coefficients in u). 



2. Antisymmetric modes. 

As in subsection [TV A| , a pattern built on m arbitrary pairs of antisymmetric modes is 
marginally stable versus a m+1 pair of wavevectors. Furthermore, since the contributions of 
the quadratically resonant symmetric and antisymmetric modes, in their respective ampli- 
tude equations have opposite signs (v\ < and v 2 > 0), pure antisymmetric mode patterns 
are also unstable versus resonantly coupled symmetric and antisymmetric ones. Recall also 
that a pattern built on m arbitrary pairs of antisymmetric modes is marginally stable versus 
non resonant symmetric modes. 



3. Mixed modes. 

As a result, when quadratic resonances between symmetric and antisymmetric modes are 
possible, pure steady patterns built on symmetric or antisymmetric modes only, are always 
unstable. We have thus to consider the possibility for the system to develop mixed modes 
patterns. 

Let us consider the simplest case of mixed mode patterns built on one symmetric mode 
and two quadratically resonant antisymmetric modes (e.g. Si, Ai and Bi in figure 0). Their 
uniform amplitude equations are: 

S, = (21. - 1)S 1 - I 8 A^ - S x [(2u 2 + u 3 )N 2 + 2w 2 (\A 1 \ 2 + \B l \ 2 )] 

A 1 = {2I S - 1)A X + hS^ - A 1 [2{u 2 - ws)^] 2 + (2w 2 \A 1 \ 2 + (2w 2 + ^i)^! 2 )] 

B x = (21, - + I S S 1 A 1 - B 1 [2(u 2 - w 3 )\S 1 \ 2 + (2w 2 \B l \ 2 + {2w 2 + w^A^ 2 )] . (26) 
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Because of the symmetry between A\ and B\ we look for solutions with the same amplitude 
for both antisymmetric modes. Defining amplitude and phase variables as S\ = R s expi<j), 
Ai = R a exp tip, Bi = R a exp i-i/j, and \I/ = <fi — ip — ip, one has 

R s = {21 s - l)R s - I s R 2 a cos * - R s [{2u 2 + u 3 )R 2 s + Aw 2 R 2 a ] 

R a = (2I S - l)R a + I s R a Rs cos # - R a [2{u 2 - w 3 )R 2 + (Aw 2 + 2w 1 )R 2 a \ 

* = ^[^-2^]sin*. (27) 

A phase stable steady state corresponds thus to \l/ = 2mr if R 2 a < 2R 2 S) and to \I/ = (2n + l)n 
if R 2 > 2R 2 . Combining the steady state conditions for R s and R a , it can be seen that 
\1> = 2rm requires that R 2 a > R 2 and that \1/ = {2n + 1)tt requires that R 2 a < R 2 . As a result, 
a stable steady state may only be obtained for \P = 2n7r, with R 2 < R 2 < 2R 2 . Furthermore, 
this condition is satisfied if the kinetic coefficients are such that 1 < |0| < 2( v / 3 - 1) ~ 1.46, 
which is thus a necessary condition to be satisfied to obtain such mixed modes solutions. 
This condition corresponds to 0.57A; S < k a < 0.65/c s or 0.64 < cos0 < 0.769, with I s ~ 0.5. 

A similar analysis may be performed for a pattern formed by an hexagonal planform of 
symmetric modes and their quadratically resonant antisymmetric ones (cf. fig0). However, 
in this case, the fact that the quadratic couplings between symmetric and antisymmetric 
modes have opposite signs does not allow the stabilization of critical patterns. 



C. Summary of the analytical results. 

The conclusion of the analysis presented is the previous two subsections is as follows: 

• For k a < 0.5k s (\6\ < 0.666), close to threshold hexagonal symmetric patterns are 
stable. For slightly larger values of I s symmetric stripes and hexagons are bistable 
and finally only the stripes remain stable. Patterns built on an arbitrary number of 
antisymmetric modes are neutrally stable. As this results come from an expansion up 
to cubic nonlinearities we can not conclude about the stability of the antisymmetric 
patterns. It would be necessary to go to higher order terms. Finally, there are no 
mixed stationary patterns. 

• For 0.5k s < k a < 0.57k s (0.666 < \8\ < 0.963), pure symmetric and antisymmetric 
mode patterns are unstable. Also, no steady mixed mode patterns are found either. 
For that reason one expects to find time dependent structures which involve both 
symmetric as well as antisymmetric modes. 

• For 0.57k s < k a < 0.65fc s (0.963 < \6\ < 1.46), one may expect steady patterns 
formed by the superposition of one symmetric mode and its quadratically resonant 
antisymmetric ones. We should notice that the range in which this steady mixed 
patterns does exist may be in fact smaller than 0.963 < \9\ < 1.46. This is because 
our calculation was for a necessary condition for the stability of the global phase, but 
this is not a sufficient condition for the stability of the pattern. 
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• For 0.65k s < k a (1.46 < \8\), symmetric modes are unstable versus antisymmetric 
modes and no steady mixed mode solutions are found. Patterns composed of an 
arbitrary number of antisymmetric modes are neutrally stable, so it would be necessary 
to include higher order nonlinearities to conclude about their stability. 

The present analysis provides the basic elements for the study of pattern formation when 
both symmetric and antisymmetric modes become simultaneously unstable. As it will be 
discussed below, it is partially confirmed by numerical analysis of the complete dynamical 
model. It could be improved, on the one hand, in determining the full stability range of mixed 
mode patterns, and, on the other hand, in resolving the issue of finding asymptotic states in 
cases where no steady critical patterns are found. In such cases, as suggested by numerical 
analysis, one should consider the possibility of asymptotic time-dependent or non critical 
patterns. The latter case would require to incorporate in the dynamics harmonics or non 
critical modes. We should notice that a peculiarity of the situation considered here is that 
there are no cubic nonlinear terms involving only critical modes. The cubic nonlinearities 
we have considered are generated solely by the adiabatic elimination of stable modes from 
the quadratic terms. However for an input field which is slightly above threshold there will 
be a range of unstable modes around the critical one. The cubic nonlinearity N 3 for these 
modes will be small but non zero. As the cubic nonlinearities generated by the adiabatic 
elimination of stable modes from the quadratic terms are also small, it would be necessary 
to include in the analysis both kinds of cubic nonlinearities. This is beyond the scope of the 
present paper. 



V. NUMERICAL RESULTS 

We have performed several numerical integrations of Eq. [IJ using a numerical scheme 



described in detail in Ref . |U| . The method is pseudospectral and second-order accurate in 
time, and is similar to the so-called two-step method. Lattices of size 128 x 128 were used. 

For \6\ < 0.666 there is not a mixed mode pattern since quadratic resonance is not 
possible due to the fact that k a < k s /2. An example of this situation is shown in fig. [| for 9 = 
—0.5. According to the weakly nonlinear analysis stationary symmetric hexagons are stable 
close to threshold, however we find numerically a disordered structure for both polarization 
components of the field E x and E y . As E x = (E + + E_)/y/2 and E y = (E + — EJ)/y/2i, the 
structure observed in E x is formed by symmetric modes with wavevectors in the ring \k\ = k s 
whereas the structure in E y is composed of antisymmetric modes with wavevectors \k\ = k a . 
The overall structure is built on an arbitrary number of symmetric and antisymmetric modes. 
The weakly nonlinear analysis predicted that there were no stationary mixed patterns and, 
in fact, the pattern shown in fig. [5] from numerical integration is not static, it evolves 
dynamically in a slow time scale. 

For 0.666 < \6\ < 0.963, although resonance is possible, no mixed mode pattern is 
numerically found, as predicted analytically. As in the previous case, a disordered structure 
is obtained and no steady pattern is reached even for long time integrations. Rings of radius 
k s and k a in the far fields of the linearly polarized components x and y respectively, show 
that there is not a selection process of critical wave vectors, all directions become unstable. 
In fig. |] we show the field configurations for 9 = —0.83. 
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The analysis of the previous section showed that mixed mode patters may be expected 
for 0.963 < \9\ < 1.46. In fact, for 9 = — 1 a pattern formed by the superposition of a 
symmetric mode and the quadratically resonant antisymmetric ones is obtained, see fig. |5|. 
A roll pattern is seen in the x linearly polarized component with wave number k s , while a 
rectangular pattern appears in the y component with wave number k a . As for the range 
of existence of mixed mode patterns, we have found numerically the existence this patterns 
also for 9 = —1.05 but not for \9\ > 1.1. 

For 9 = —1.1 we find a steady pattern in which the symmetric modes are damped 
and a stripe pattern with wave number k a appears in the ^-polarized component, while 
x-polarized component is almost homogeneous (in fact, a small amplitude stripe pattern of 
wavevector 2k a is observed in E x which comes from the coupling between E x and E y ). This 
pattern made out of antisymmetric modes is the same pattern that is formed when only 
the antisymmetric instability is present. Overall, the electric field displays an elliptically 
polarized spatial structure. Similar stationary patterns are found for larger absolute values 
of the detuning. In fig. |6| we show the case with 9 = —2; this value was chosen in order 
to induce the formation of a square pattern, since the angle <ft between symmetric and 
antisymmetric vectors is <fi = 45° in this case, but anyway the antisymmetric stripe pattern 
is formed as described. 

VI. CONCLUDING REMARKS 

We have studied the spatial polarization structures in a mean field model for a Kerr 
medium close to a two-photon resonance and driven by a linearly x-polarized field with 
negative detuning. In the self-focusing case the first pattern that is formed as the pump 
intensity is increased arises from the competition between two stationary instabilities which 
occur simultaneously. This codimension two bifurcation appears naturally in the system and 
it is a consequence of the form of the nonlinearities associated to the two-photon-resonant 
four-wave mixing process. As this codimension two bifurcation is not the result of the 
fine tuning of two parameters, as it is usually the case, it should be much more simple 
to be observed. In fact there is only one parameter, the pump, which has to be tuned. 
Furthermore, we still have another free parameter, the detuning, which allows the system 
to form different patters without changing the distance to the codimension two instability 
threshold. 

Near the instability threshold, we have obtained the evolution equation for the patterns 
arising from the interacting instabilities using a weakly nonlinear analysis. From these anal- 
ysis and from the numerical integration of the model we have shown that we can have the 
following patterns: a) dynamical structures involving an arbitrary number of symmetric 
and antisymmetric modes, b) a steady state pattern formed by the superposition of a sym- 
metric mode and the quadratically resonant antisymmetric ones and c) an antisymmetric 
stripe steady pattern in E y which is the same that would appear if only the antisymmetric 
instability were present. 

As a final remark we notice that an interesting peculiarity of the situation considered 
here is that the selection among these patterns can be done in a very natural way, that is, 
changing the value of the detuning. 
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FIGURES 
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FIG. 1. Marginal stability curves for linearly polarized input field corresponding to the sym- 
metric solution. In the inset the symmetric steady state homogeneous solution is shown, as a 
function of the input field intensity, for linearly polarized light. Value of the detuning 9 = —1. The 
parameters used in Eq. (|l]) in order to get the codimension two situation are: r) = l, a = l, A = 
and B = 2. The quantities plotted in all the figures are dimensionless. 
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FIG. 3. Field configuration for 9 = —0.5 and I s = 0.51 (Jo = 3.31) after integrating Eq. |l] for 
a time t = 35000. From left to right and from top to bottom: $t(E x (r)) (near field x-polarized 
component plotted with the grayscale: black=0.51, white=0.60), \E x (k)\ 2 (far field x component), 
?R.(Ey(r)) (near field y-polarized component plotted with the grayscale: black=-0.035, white=0.042) 
and \E y (k)\ 2 (far field y component). 
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FIG. 4. Field configuration for 9 = —0.83, I s = 0.51 (Jo = 4.41) and an integration time 
of t = 65000. From left to right and from top to bottom: $l(E x (r)) (grayscale: black=0.44, 
white=0.51), \E x (k)\ 2 , %t(E y (r)) (grayscale: black=-0.036, white=0.033) and \E y (k)\ 2 . 
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FIG. 5. Field configuration for 9 = — 1, I s = 0.51 (Jo = 5.06) and an integration time 
of t = 14000. From left to right and from top to bottom: $l(E x (r)) (grayscale: black=0.43, 
white=0.47), \E x (k)\ 2 , ^k{E y {f)) (grayscale: black=-0.047, white=0.047) and \E y (k)\ 2 . 
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FIG. 6. Field configuration for 9 = —2, I s = 0.51 (Jo = 10.2) and an integration time 
of t = 10000. From left to right and from top to bottom: $l(E x (r)) (grayscale: black=0.32, 
white=0.34), \E x (k)\ 2 , ^{E y {f)) (grayscale: black=-0.15, white=0.15) and \E y (k)\ 2 . 
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